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We show the existence of a very large spin Hall effect of light (SHEL) in single photon plasmonics 
based on spontaneous emission and the dipole-dipole interaction initiated energy transfer (FRET) on 
plasmonic platforms. The spin orbit coupling inherent in Maxwell equations is seen in the conversion 
of (T+ photon to (j~ photon. The FRET is mediated by the resonant surface plasmons and hence we 
find very large SHEL. We present explicit results for SHEL on both graphene and metal films. We 
also study how the splitting of the surface plasmon on a metal film affects the SHEL. In contrast to 
most other works which deal with SHEL as correction to the paraxial results, we consider SHEL in 
the near field of dipoles which are far from paraxial. 



In recent years the SHEL has attracted considerable 
attention \^\^. Several experiments have confirmed the 
existence of the effect [SHlOj. The simples t versio n of 
the effect is known as Federov-Imbert effect 
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is seen prominently as a polarization dependent trans- 
verse shift on the reflected beam at a dielectric interface. 
The shift occurs relative to the prediction of geometrical 
optics. The light beam is known to carry both orbital 
and spin angular momentum [l4, 15 1 and the spin-orbit 
coupling is inherent in the vectorial Maxwell equations. 
When light travels across the interface between the two 
dielectrics then the orbital angular momentum changes 
which then implies change in the polarization so that the 
total angular momentum is conserved [l6l |. The mag- 
nitude of the shift is generally quite small. However, 
for reflection from metal surfaces shifts of the order of 
wavelength have been reported Q. The metallic grat- 
ings even yield larger result [1] due to excitation of sur- 
face plasmons. The SHEL appears in a variety of other 
ways iig. 

Motivated by the recent progress in the observation of 
very large enhancement in spontaneous emission [l^1 - 124] | , 
we propose the existence of giant SHEL in single photon 
plasmonics. Especially fabricated plasmonic structures 
(PS) [2^ 26 1 have been shown to be especially suitable 
for enhancement of fluorescence. We consider the Forster 
energy transfer 27, 2^ between two atoms located on a 



plasmonic surface. The energy transfer is mediated by 
plasmons. The SHEL arises as we demonstrate as a clear 
conversion of a, say, (t"*" photon into a a~ photon which 
is absorbed by the second atom leading to preparation of 
the second atom in a state which is orthogonally polar- 
ized to the state of the atom which produced the photon 
in the first place. The conversion is accompanied by the 
generation of the two units of orbital angular momen- 
tum. The giant SHEL can be detected by monitoring the 
population of the orthogonally polarized excited state of 
the second atom. We trace the existence of the effect 
to the fact that the dipole field is far from a paraxial 



field and in fact contains all the Fourier components in- 
cluding the evanescent waves. It has been shown that the 
image of a dipole in far field is displaced 13 , 2^ . Further 
Maxwell equations imply that the vectorial properties of 
the fields change on propagation. We demonstrate giant 
SHEL with explicit results for FRET on graphene and 
metallic platforms. The atoms on nano fibers [13, Ej 
would be another prominent system which would lead to 
large SHEL. Our proposal constitutes a new generation 
of possible experiments for large SHEL. 

Let us consider the radiative interaction and the en- 
ergy transfer between two three-level atoms A and B 
as depicted in Fig. [TJ The transfer of the excitation of 
atom A, say, to atom B occurs via spontaneous emission 
from atom A and absorption of the emitted photon by 
atom B. Thus the single photon excitation transfer is 
due to the interaction with the vacuum field E in pres- 
ence of the plasmonic environment. For concreteness we 
assume that atom A has been prepared by pulse excita- 
tion in the excited state \ep^,m = +1). Then atom A will 
emit a polarized photon and drop to the ground state 
I (7a)- Clearly this (t+ photon can excite the transition 
I (7b) — \eB,,rn = +1). We use a fixed coordinate system 
as shown in Fig. [1] The quantization axis is taken to be 
perpendicular to the PS. We introduce the circular polar- 
ization vector e± = (e^jiiej,)/ V2 where and ey are the 
unit vectors in x and y direction, respectively. Now, we 
ask the question — can the photon emitted by atom 
A excite the transition jge) ^ Icb, = ^1) in atom _B? 
Clearly this requires a a~ photon defined with respect 
to the fixed coordinate system. This would be possible if 
the propagation of the electromagnetic field is such that 
the photon gets partially converted into a cr^ photon 
when it reaches atom B. This conversion would be due to 
the spin-orbit interaction inherent in Maxwell equations. 
Let us denote by I the orbital angular momentum then if 
is converted into a~ such a conversion would happen 
by producing two units of I, i.e. |cr+, Z = 0) — > ,1 = 2). 
Thus the excitation of atom B to the state \e^,m = —1) 
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is the SHEL in FRET. It should be borne in mind that 
conversion of to a~ would happen with a certain prob- 
ability and hence the probability of exciting atom B to 
the state \e^,m = —1) is nonzero. The SHEL in FRET 
can be monitored by studying the populations of the ex- 
cited states by applying another field and by studying far 
off resonant flourescence. 

We next present our method of calculation. The 
Hamiltonian in the interaction picture can be written as 
in the interaction picture by the Hamiltonian 

Hi{t) = -^p„j •E(r„,t)|e„,,)(g„|c'"«* +h.c. (1) 

Here the sum is over both the atoms (a = A, B) and 
the two excited states of each atom. Note that (e, m = 
±l|p|5) = pe*^ where p is the amplitude of the dipole 
matrix element. Further E is the Hermitian electro- 
magnetic field operator. We have four possibilities of 
energy transfer: two with Am = (|eA,TO = +1) 
\eB,m = -1-1), \eA,m = -1) \eB,Tn = -1)) and 
two with Am = ±2 (|e^,m = —1) |es,m = +1), 
\eA,m = +1) \eB,m = -1)). 

If we assume that atom A is prepared in an excited 
state as depicted in Fig. [T] (a) then the wave function of 
the system at time t is up to second order in Hi given by 



dh / dt2 Hi{h)Hiit2)\ibo) (2) 



m)) = 



where |V'o) = |e^, m = -1-1, .gs, {0}); {0} denotes the vac- 
uum state of the electromagnetic field in the presence of 
the PS. 

We are now interested in the probability of exciting 
atom i? at a later time. Therefore we need to calculate 
the transition amp litude C ~ ((^a, e^, {O}|'0(i))- It can 
be shown that [321 



dC_ 
dt 



• G(rB,rA,a;o) • (3) 



where the asterisk symbolizes the complex conjugation. 
Here, we have introduced the dyadic Greens function 
G of the environment. Note that the full calculation 
leading to Eq. ([3]) includes both resonant and nonres- 
onant quantum paths |eA,5B) |5a,.9b) |.9A,eB); 
|eA,5B) |eA,eB) -> |gA,eB). Eq. © gives the prob- 
ability which can be converted in the usual manner to 
the transition rate by summing over the final states. We 
normalize the transition rate by the Einstein A coefficient 

r = for the level |eA,m ~ +1). We call this ratio 

as Dn [for FRET to |eB,m = +1)] and D^i, [for FRET 
to [cb, m = -1)] 
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(Gj;^; — Gyy + i[G2-y + Gj^j;]) 
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(4) 
(5) 



It is important to note that we have finite probabilities 
for both kinds of excitation transitions which refer to a 
zero spin-angular momentum coupling (Am = 0) or a 
non-zero spin-angular momentum coupling (A?ti ~ ±2). 
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Figure 1: Sketch of the excitation transfer between atom A 
and B. 

Let us first consider the process of FRET in free space. 
Then we can determine the transition probabilities by 
means of the free Green's function G^"-* which is given by 
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(6) 



When assuming that both atoms are in a plane perpen- 
dicular to the z axis defined by the quantization axis and 
let d be the interatomic distance, then we obtain 
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■6e+2'^ , (7) 
(8) 

where we have introduced a = [k'^d^ + ifcorf — \)/k'^d'^ , 
and b ~ {?> — iik^d — k'^d^) / k^d^ \ fco ~ liJqIc- The angle tp 
appearing in the first expression is the angle between the 
line connecting the positions of the atoms A and B with 
the X-axis. Eq. ([7]) shows that during the excitation tran- 
sition -Dgjj the orbital angular momentum has changed 
by two units which is a manifestation of the spin-orbit 
interaction. 

To see which of both excitation transitions is more 
likely, we plot in Fig. [2] the relative transition probability 
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for free space. It can be seen that for the considered 
wavelength region and the interatom distances d = 30 nm 
and d = 50 nm the SHEL is prominent, i.e. energy trans- 
fer with change in spin-angular momentum is generally 
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speaking more likely than the one with no change in spin- 
angular momentum. Note also that for small distances 
kod <C 1 we have a w — l/(fcod)^ and b w 3/(fcod)^ so 
that i? 9; for large distances kod ^ 1 we find i? w 1. 




1000 



Figure 2: R from Eq. ([9| as function of the wavelength. The 
inset shows a semi-logarithmic plot of from Eq. ((TJ for 
the same wavelength range. 

Now, we want to show that the SHEL in free space can 
be enhanced up to several orders of magnitude by the 
interaction of the two atoms A and B with the surface 
plasmons of a given sample. To this end, we first need to 
determine Dgh and -D„ for the case where the atoms are 
placed in a given distance z above a PS (see Fig.[T]). This 
can be done by using a spacer material on the sample, 
for instance. Then the Green's function is given by the 
sum G = G^'"'^ + G'^'*^) of the vacuum Green's function 
G'^"^ and the scattered Greens function G^'^'^'' which takes 
the interaction with the surface into account. Analytical 
expressions for G'-*'^-' for a sample with a plane surface 
can be found in Ref . [s^l ■ When inserting the full Greens 
function into Eqs. (U) and ([S]) we obtain 
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tions. and rp are the amplitude reflection coefficients 
of the PS for s and p polarization. The poles of rp de- 



termine the surface modes of the PS. Further, Eqs. (|T0|) 
and (jlip are far from paraxial, since all possible k are 
taken into account: the propagating waves {n € [0, fco]) 



as well as the evanescent waves (k G [fco, oo[). This time 
the coupling to the angular momentum manifests itself 
not only through the phase factor e^^"^ but also through 
the different Bessel functions. 

Finally, we introduce the plasmonic enhancement fac- 
tor 



E = 
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which allows us to quantify the enhancement of the SHEL 
by means of the interaction with the surface plasmons of 
the sample. As a first example let us consider a sus- 
pended sheet of graphene. The reflection coefficient for 
the p-polarized modes is in this case given by [s^] 



47rf77Q 
, 47rg7o 



(13) 



Using this quantity with the analytical expressions for 
the conductivity cr(w) of graphene [34| we obtain the plas- 
monic enhancement E plotted in Fig. |31 In Fig. EJa) one 
can observe the frequency dependence showing a broad- 
band effect of plasmonic enhancement due to the interac- 
tions of the plasmons in the graphene sheet which exist 
for frequencies smaller than 1 eV for the chosen doping 
level. In Fig. ^h) we show the distance dependence of 
the plasmonic enhancement choosing different frequen- 
cies. The curves are normalized to the corresponding 
propagation length of the surface plasmons I = l/Im(Ksp) 
which can be derived from the poles of the reflection co- 
efficient in Eq. ([T3| . Obviously, the enhancement dies 
off for interatom distances d ^ I which shows that the 
plasmonic enhancement is due to a excitation transition 
mediated by the surface plasmons. Note, that for certain 
distances the plasmonic enhancement leads to a transi- 
tion probability which is six orders of magnitude larger 
than for vacuum. That the SHEL is so large can be ex- 
plained by the fact that due to the interaction with the 
surface plasmons of the PS the emission of the photon by 
atom A is enhanced as well as the absorption by atom 
B. Thus local field enhancement factors contribute twice 
making th SHEL giant. 

Next we consider the role of surface plasmons on a 
semi-infinite metal surface. In this case the reflection 
coefficient is given by the usual Fresnel coefficient 



7oe - 7 
7oe + 7 



(14) 



where we have introduced 70 = yfcj— the wavevec- where e 
tor in z direction, k = (kxykyY the wavevector parallel 
to the surface; Jq and J2 are the cylindrical Bessel func- 



the permittivity of the sample and 7 = 
^/fcge — K^. Here, we use a Drude model 1351 which has 
been fitted to the Johnson-Christy data |36|- The re- 



sulting plasmonic enhancement E is plotted in Fig. |4l 
The wavelength dependence in Fig. IH^a) shows a large 
enhancement for wavelenght of Xsp = 291, 7nm. This is 
exactly the surface plasmon resonance frequency of silver 
which follows from Rc[e(a;sp)] = — 1 and Im[e(ajsp)] ^ 1. 
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Figure 3; (a) The giant SHEL on suspended grapliene with 
Ef = 1 eV as function of wavelength, (b) E as function of 
distance. The distance is normalized by the corresponding 
propagation length of the surface mode / — 1/Im(fi:sp) which is 
1(0.2 eV) = 2138nm,Z(0.3eV) = 1392 nm, Z(0.5 eV) = 772 nm, 
and 1(0.7 eV) = 486 nm. 



Here, the observed plasmonic enhancement results in a 
probability Dgh being four orders of magnitude larger 
than for vacuum. In addition, this is a narrow-band ef- 
fect compared with graphene which is due to the fact 
that for the chosen distance z we are in the electrostatic 
regime of the surface mode. We have also determined the 
pure pole contribution to D^h which describes the narrow 
peak in Fig. 4(a). Note also the Fano-like line shapes for 
SHEL which are due to the interference of the contribu- 
tions from the pole of Vp and the non-pole contribution in 
Eq. pO|) . In Fig.m^b) we plot the plasmonic enhancement 
factor as a function of the interatom distance d. Again, 
each curve is normalized to the propagation length of 
the surface plasmons in silver I = l/Im(Ksp) which is 
now obtained from the poles of Eq. p^ . Once more, one 
can observe that the plasmonic enhancement dies off for 
d > L 

As in studies of spontaneous emission 
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the possibility of tayloring the properties of the sur- 
face modes will affect the plasmonic enhancement of 
the SHEL as well. This can for example be done by 
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Figure 4: (a) E for semi-infinite silver as function of wave- 
length. The inset shows Dsh/D„ for the same distances and 
wavelength range, (b) E for semi-infinite silver for different 
wavelengths close to the surface mode resonance as function 
of distance d. The distance is normalized the corresponding 
propagation length of the surface mode I = 1/Im(fi:sp) which 
is /(300,nm) = 692 nm, Z(305, nm) = 1458 nm, Z(310,nm) = 
2436 nm, and Z(315,nm) = 3610 nm. 



sputtering thin metallic films on a dielectric substrate. 
Then the surface mode branch will split into two sur- 
face mode branches [s^l- This splitting should also be 
obsverable in the Spin-Hall effect. In Fig. [S] we show 
the plasmonic enhancement factor for two atoms above 
a substrate made of a 15 nm silver film on a silica sub- 
strate (fsubstratc — 2.1025). The splitting of the surface 
plasmons can be nicely seen. Furthermore, it can be 
observed that the excitation transfer mediated by the 
long range plasmon (s^l (A < 300 nm) leads to a much 
larger enhancement than that by the short range plasmon 
(A > 320 nm) which is similar to the enhanced emission 
rate of a molecule close to a corrugated metal film [l^ . 

In conclusion we have presented a new way to real- 
ize the spin Hall effect of light which is inherent in the 
near field of an oscillating dipole. We showed how this 
kind of SHEL can be monitored by observing the energy 
transfer to an orthogonally polarized atomic state. The 
change in polarization is accompanied by the generation 
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Figure 5: E for a 15 nm silver film on silica as function of 
wavelength. The inset shows a 3-d plot of Iog(lm[rp(ctj, k)]) 
from Eq. (|14|l for the same wavelength range and lateral 
wavevectors ti/ko £ [1, 10]. 



of two units of orbital angular momentum. The SHEL 
is largely enhanced on plasmonic platforms due to reso- 
nant surface plasmons. Other systems like nano fibers or 
nano antennas could be equally useful in producing large 
SHEL. 

The idea was initially presented by GSA at the TIER 
school on Plasmonics 2012 in Hyderabad and GSA 
thanks the Director TIER for hospitality. 
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